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ABSTRACT 

An inert sphere of a few meters diameter, placed in a special stable geosynchronous 
orbit in perpetuo, can be used for a variety of scientific experiments. Ground-based 
observations of such a sphere, “GeoSphere”, can resolve very difficult problems in 
measuring the long-term solar irradiance. GeoSphere measurements will also help us 
understand the evolution of Earth’s albedo and climate over at least the next century. 
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1 MOTIVATION 


The Sun is a very stable object. The Sun’s “irradiance” 
(bolometric flux measured since 1978 from above the earth’s 
atmosphere) varies by 0.06-0.1% peak-to-peak, on time 
scales of a 11 years (see, fo r example, recent reviews by 
iFrohlicl] l2013l : IWillsod 1201^. This vari ation follows the 
well known “sunspot cycle” ( Schwabell844l l. On shorter time 
scales, the evolution and disk passage of sunspot groups pro¬ 
duce variatio ns that are often 0.2%. Given that sunspots 
are magnetic llHaldflOO^ l. these variations arise because of 
a variable solar magnetic field. While of small amplitude, 
these changes are far more rapid than expected from first 
principles. They are believed to be generated by a magneto¬ 
hydrodynamic dynamo, drawing energy from differential ro¬ 
tation and turbulence in the Sun’s interior. The dynamo 
engine changes the Sun from a benign object, stable over 
time scales of 100,000 years and radiating essentially as a 
black body near 6000K, to a rapidly varying object.The dy¬ 
namo also c auses far mor e variable high energy UV and X- 
ray photons (IWhitdll977l l. wind variations, and intermittent 
energetic particles and magnetized plasma clouds. 

The earth’s climate is sensitive to several factors, the 
Sun is by far the largest external factor. Herein lies a prob¬ 
lem. We have no physical explanation of the solar dynamo, 
thus our understanding of irradiance changes remains em¬ 
pirical. Nor do we have accurate measurements of solar vari¬ 
ability on time scales longer than decades. Long-term vari¬ 
ability is a contentious and active research area, plagued by 


the absence of hard data fe.g. ljudge et al.ll2012lL Irradiances 
have been measur ed from 1978 using radiometers in space 
le.g. IWillsonlbOl^l . Pre-flight calibration precisions of «500 
parts per million (ppm), 0.05%, have been achieved. But in¬ 
dependent data from different space missions have failed to 
provide the needed long-term stability. First, the highly pre¬ 
cise irradiances measur ed by differen t experiments differ by 
0.5-1% (see figure 1 of IWills^l2014ll . i.e. disagreements of 
10 to 20 (t, on timescales longer than 11 years. These system¬ 
atic thwart attempts to measure variations longer than one 
satellite lifetime. Second, instrumental sensitivities change 
over mission lifetimes which are, unfortunately, rather close 
to 11 years. Third, expensive new radiometers must be reg¬ 
ularly launched. 

We propose an alternative solution that can provide 
society with solar brightness measurements that are stable 
over centuries at a time when climate change is of first im¬ 
portance. Our long-term precision photometric goal is set 
to 0.01% averaged over a week. This may not be achiev¬ 
able, thus our less stringent requirement is a 0.1% precision 
over a century, which will uniquely constrain long-term so¬ 
lar va ri ations at a l e vel st ill of great interest IShapiro et al.l 
ll201lll ; iLubin et al.l (l2012l ~) . We adopt high precision differ¬ 
ential photometry of a target in space against stable en¬ 
sembles of stars. Nightly precisions of stellar observations 
of “only” 0.1% are achieved w ith existing small telescopes 
l|Young et al.lll99ll :[H enr vI 1999D . with seasonal averages pre¬ 
cise to 0.01%. 
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2 A SIMPLE SOLUTION 

Photometry of planets, asteroids and the moon present dif¬ 
ficulties when trying to extract solar variations near 0.01%. 
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Planets have atmosphere-rela ted albedo variations ( for Nep¬ 
tune, mv ~ 8 , they are « 4%. lLockwood et al.lll99ih . Aster¬ 
oids have variable brightnesses dependent on their rotational 
aspect and highly variable geocentric distances. The moon 
is so bright and extended, with diverse phase-dependent 
brightnesses of its features, that accurate relative photome¬ 
try with standard stars is not possible. 

We suggest that a diffuse object, “GeoSphere”, be 
placed in a stable geosynchronous orbit (GSO). The need 
for a high level of geometric symmetry argues strongly for 
an unemcumbered sphere. A sphere of a few meters diam¬ 
eter has mv ~ 10 at opposition independent of orientation, 
and the optical appearance of a star (Table 1. Placed at 
L 2 , mv ~ 18). Photometric observations of GeoSphere with 
ground-based telescopes can in principle provide century- 
long monitoring of GeoSphere’s brig htness at 0.01% preci¬ 
sion dYoung et al.lll99ll : |h enr using statistical aver¬ 

ages of many (« 20-I-), nightly observations. If the scattering 
properties of the sphere in orbit can be measured, calibrated 
and/or modeled at the same level of precision, solar bright¬ 
ness variations will follow. GeoSphere should spin with a 
period of a few seconds, less than a typical integration time 
of precision photometry (a minute). 

Table 1 lists conditions for a 4m sphere as seen from 
Earth. We adopt a 4m diameter only because fairings of 
large launch vehicles extend to 5m. Smaller vehicles might 
be used with on-orbit expansion of a pre-com pacted spher e 
(e.g., a “Hoberman sphere”). Using data from I All^ (Il973ll . 
the Johnson mv-magnitude when the sphere is at angle a 
from the anti-Sun point, is, with radius Vs = ds/2: 

mvia) « -26.74 - 2.5logio (7 - ——■ (1) 

\{ro-r(syq J 

Here we use a perfectly diffusing sphere whose phase func¬ 
tion 4>{a) obeys Lambert’s law: 

(/(q) = [sin a-\- [tt — a) cos a)] /tt (2 ) 

for which q = 1.5. Observing GeoSphere close to opposition 
give s (j>{a 1 ) ~ 1 , and mv = 10.08. 

lYoung et al.l lIlOOll . their table 1) calculate exposure 
times needed to achieve 0.05% RMS photon noise in the 
Stromgren b and y filters at 4670 and 5470 A, with widths 
of 160 and 240 A respectively. Assuming mv [GeoSphere) = 
10.08, an integration time of < 100 s, we find that telescope 
aperture D and GeoSphere diameter ds must satisfy: 

[dsD)^ > 36 m"^ (Stromgren by at opposition). (3) 

In the general case of a spectral sampling of w A, 

[dsD)^ > 5800/u; mA (4) 

For observations with the w = lOA s ampling of the SORCE- 
SIM instrument llHarder et ahllioo^ j for example, then ds = 
4m and D > 6 m. 


3 SOLAR BRIGHTNESS OBSERVATIONS 

Narrow-band photometric measurements of solar-like stars 
are routinely made using small (< 0.8 ) Automated Ph oto- 
electric Telescopes lAPTs. [Henr ^19991: iHall et al. I I 2 OO 9 II . In 
practice, combining all sources of error, those brighter than 


Table 1. Conditions for a 4m Geostationary sphere 


Parameter 

value 

Environmental conditions at GeoSphere 

Geocentric distance Vo km 

42,164 

Solar my 

-26.74 

Earth my (/3 = 90°,quadrature) 

-20.6 

Earth my (/3 = 20°) 

-16.0 

Moon my 

> -12.74 

Physical and observable GeoSphere properties 

Adopted sphere diameter ds m 

4 

Angular diameter arcsec 

0.0196 

Adopted albedo a 

0.9 

q (phase law parameter) 

1.5 

my (opposition) 

10.08 

my (quadrature) 

11.32 


Notes: “Lambert’s law” for scattered radiation is assumed, corre¬ 
sponding to a perfectly diffusing sphere. Magnitudes were com¬ 
puted assuming the distance to the sphere is To— r^. is the angle 
between the geocenter-GeoSphereline and geocenter-Sun line. 


mv = 8 have a daily rms precision of 1 milli-magnitude in 
b and y. Scaling to the brightness of GeoSphere, telescopes 
of 1.5 m diameter are needed to measure the brightness of a 
4m diameter GeoSphere at the same daily level of precision. 
With telescope time devoted exclusively to GeoSphere, the 
precision will improve over the 0 . 1 % nightly means reported 
in the stellar APT programs. 

However, qualitatively different observations must be 
made to permit differential photometry of GeoSphere 
against stars. Objects in GSO are observed most simply 
by switching off a telescope’s astronomical drive. In time 
t the Earth’s rotation sweeps out 15” cos 5 t of sky, where 
S is the object’s declination. As a result, background stars 
passing near to GeoSphere will contribute to the integrated 
flux. Gonsider a photometer fed from an aperture in the 
telescope’s focal plane of 5” diameter, 5 = 0°, and t — 100 
seconds. The area covered is 6 x lO”"^ square degrees. If 
the average number of stars with visual magnitude > m is 
N[m,b) per square degree at Galactic longitude b, then the 
number of stars crossing the field per lOOs observation is 

n = 6 X 10~'^N[mv,b)- 

Using data from (lAllenI Il973l . §117), for mv = 10 , 

N{mv,b) varies between 1.25 (Galactic latitude b = 
0 °) and 0.55 [b = 90°); for mv = 20 we find logjQ N[mv, b — 
0,90) = 5.0 and 3.4 respectively. For 6 = 0 (worst case), 60 
stars of magnitude 20 or brighter will cross the aperture; for 
mv < 10 , on average 1 in 100 such integrations will have 
one star of brightness similar to GeoSphere affecting the 
integrations. One solution might be to feed three mutually 
calibrated photometers from the focal plane of the telescope, 
offset along a line of constant declination. With GeoSphere 
measured from one part of the focal plane, the same stars 
can be measured within a second before and after the tar¬ 
get in the other photometers, allowing subtraction of their 
contribution. 
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4 EARTH ALBEDO VARIATIONS 

Changes in Earth’s albedo might be monitored observing 
GeoSphere close to quadrature, where it is illuminated by 
the Earth’s sunlit hemisphere. Then ao < tt/ 2, 4>{oio) 
I/tt and the brightness of GeoSphere from dire ct sunlight i s 
mv > 10.96. Using mean earth albedo data (|Allenl[l97,lh . 
scattered light from the Earth adds « 0.4% of the direct 
solar flux to the GeoSphere at quadrature. By measuring 
GeoSphere’s differential brightness as it crosses the sky, with 
statistical averaging, the goal of of 0.01% precision might 
be achieved. Long-term relative changes in Earth’s albedo 
greater than 2.5% should be measurable by GeoSphere. 


5 TECHNICAL CHALLENGES 

5.1 Orbital stability 

The orbits of GSO satellites are perturbed by non-spherical 
components of the Earth’s gravitational field, by lunar influ¬ 
ences and by momentum transfer from solar radiation. Even 
so, certain orbits are stable enough over decades, for the pur¬ 
poses discussed here. Orbital drift through the Poynting- 
Robertson effect is unimportant for conceivable area-to- 
mass ratios for a s phere. Quoting from the conclusions of 
iFriesen et al.l (| 19921 '!: 

“Orbits of large satellites are radially stable throughout the re¬ 
gion within 2000 km above and below the geosynchronous radius. 
Satellites placed in such orbits can be expected to go no more 
than 50 km above their original apogees or more than 50 km 
below their perigees for at least several centuries, barring colli¬ 
sions.... ‘A stable plane’ exists at the GSO distance., for which 
orbit plane motion is quite limited... if the initial angular momen¬ 
tum vector of a circular GSO is aligned with the 7.3° displaced 
axis about which precession takes place for the initially equato¬ 
rial case, orbit plane motion will be limited to 1.2° of the initial 
plane.” 

A 50 km change in orbital radius gives a 1 mmag change 
in brightness. These theoretical considerations are borne out 
in part by measurements of the orbits of abandoned satel¬ 
lites. INTELSAT 1, abandoned on December 2 1969, has re¬ 
mained in an orbit with semi-diameter between 42,130 and 
42,200 km with ecc entricity between 0 and 1.6xl0~® (e.g. 
lUlivieri et al.l l2013l f for over 41 years. There are two lon¬ 
gitudes of orbital stability arising from higher order terms 
in the dynamic equations (IFriesen et al.lflO^ 'l. 105W and 
75E, running through western North America and the mid¬ 
dle of Asia respectively. These would be the desired orbits 
for GeoSphere. 


5.2 Collision hazard 


IFriesen et al.l (1 19921 ) compute the “collision hazard” for or¬ 
bits of small inclination (closing speed of two objects in GSO 
being UcSini with Vc the orbital speed of 3 km s“^): 


CH 


VcCr 

27r2R2AR 


-1 


s 


( 5 ) 


where a ~ 7r(cis/2)^ is the geometric cross section of the 
colliding objects, R is the GSO radius and AR the rms alti¬ 
tude range of the object, ~ 50 km. Per impacting object this 
becomes 1 impact in 10® years. There are only of order 400 


satellites in GSO around earth, and < 300 other small de- 
tected objects large enough for concern dKessler et ahllioill . 
Figure 2.2). The largest threat to GeoSphere is from smaller 
natural micrometeoroids. 

Gonversely, a non-powered GeoSphere represents a col¬ 
lision hazard to other spacecraft in GSO. The International 
Telecommunication Union (ITU) has mandated an end-of- 
life strategy for GSO objects, with ability to move out of 
the GSO region a requirement for new satellites. The sci¬ 
ence goals of GeoSphere are to remain in orbit for perpe¬ 
tuity, with no pre-planned en d-of-life anticipated. The “sta¬ 
ble plane” solutions found by IFriesen et al.l (Il992l ) minimize 
relative speeds between objects in these orbits, and reduce 
space debris by reducing both impact rate and debris scat¬ 
ter. If these orbits are adopted as “graveyard orbits”, there 
will be no ITU problem. 


5.3 Secular changes in scattering properties 

The GSO space environment is not “clean”. Surfaces ex¬ 
posed to propellant, silicone, and/or exposed to the solar 
UV fl ux and ionizing rad iation, change their reflective prop¬ 
erties dPever et al.ll20li ). Oxidation is a serious issue in low 
earth orbit (LEO) where oxygen densities are 10®“® cm“® in 
typical models, but oxygen concentrations in GSO are some 
10® times smaller. Instead, in GSO, high energy particles 
and photons are the larger concern. For example, the ab- 
sorptance (1 minus albedo) of white paints commonly used 
on spacecraft increases by factors of 1.3 to 2.8 after « 100 
days exposure to solar- like UV and energetic electrons in 
LEO dPever et al ] |20ll Table 4). To achieve an albedo sta¬ 
bility over a century, GeoSphere’s surface should be a sim¬ 
ple metal, alloy, or crystalline material not susceptible to 
such changes, like telescope optics in space. Such substances 
would need testing before launch, and perhaps “pre-aged” 
by deep exposure to energetic ions, electrons and photons. 

All objects in GSO develop craters, punctures and 
chemical contamination from micrometeoroid (MM) impact. 
Man-made orbita l debris is less impor tant in GSO orbits and 
is neglected here dKessler et al.l[201ll ). We consider particles 
mostly in excess of 1 micron size, those that are smaller be¬ 
ing expected to exhibit unbound orbits because of the strong 
solar radiation pressure. 

As damage accumulates the phase function of Geo¬ 
Sphere will evolve from, for example, Lambertian to 
cratered, or lunar-like. Diverse datasets collected from ex¬ 
tended space missions indicate the effects of MM impacts 
on surfaces. The large solar cells on space station MIR in¬ 
dicated an accumulated surface coverage o f MM damage 
of 0.0 45% in 10 years of exposure in LEO dSmirnov et al.l 
I2OOOI) . Measurements in GSO are rarer. But the gravita¬ 
tional focusing for typical MMs is modest (v ~ 20 km s“^) 
compared with Earth escape speed (vesc = 11 — 4 km s“^ 
in LEO and GSO respectively) and so fluxes will be similar. 
GeoSphere will accumulate MM surface damage that cov¬ 
ers « 0.45% per century, a fractional rate k = 4.5 x 10“® 
y“^. Like the moon, cratering enhances 180° scattering at 
the expense of scattering at smaller angles. Ghanges to the 
geometric albedo may be smaller, depending on the surface 
material used and the surface optical contamination. Spill¬ 
over of impact debris beyond the craters’ areas should be 
negligible if electrical potentials of attraction (polarization 
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or net charge) are less than the kinetic energy of any charged 
exploding debris ejected outwards from the sphere. 

Damage of 0.45% in area translates to smaller than 
0.45% changes in scattered light. The albedo after time t 
in years is 

at = (1 — kt)ao + {kt)b = ao — kt{ao — b) (6) 

where ao is the initial albedo and b the albedo of damaged 
areas. Since 6 > 0 MM-induced degradation in albedo occurs 
slower than (1 — kt)ao. Nevertheless, it should be mitigated, 
modeled and/or measured. One solution might be to expose 
GeoSphere to cosmic-like dust particles to mimic a suffi¬ 
ciently long exposure (> k~^ y) in geosynchronous orbit, 
such that the phase function remains constant. A hyperve¬ 
locity gur0 might be used. Another might be to calibrate 
the scattering properties of GeoSphere using future ground- 
based stable laser sources. 


5.4 Comparison Stars 

To achieve long term precisions of fractions of milli- 
magnitudes, a GeoSphere observing program would have to 
devote effort to identifying suitabl y non-varying compari¬ 
son s tars. This is not a trivial task (lHenrvlll999l : lH^all et al.l 
I2OO9II . Geosynchronous objects circumnavigate the celestial 
equator once per day, the right ascension at opposition drift¬ 
ing about 1 degree/day. GeoSphere would therefore require 
a significant number of cross-calibrated comparison stars in 
order to place its differential magnitude on a uniform scale 
throughout the year. Typical separations within photomet¬ 
ric standard groups of three stars are near 3° and always 
less than 13° ( Hall et ^l2007l : iLockwood et al.ll2007l '). As a 
rough estimate, for a maximum intra-group separation of 4°, 
we find that 345 highly stable stars distributed along the ce¬ 
lestial equator would need to be identified, with magnitudes 
nea r that o f GeoS phere. 

iHenrvI lll999ll found that stars of spectral class F0-F3 
and A8-A9 have rms probabilities of long-term (seasonal) 
stability better than 0.0005 mag of 63.3% and 87.5% re¬ 
spectively, and that short-term photometric stability tends 
to correlate with long-term stability. These findings should 
be used in any comparison star search. The observing pro¬ 
gram of the current Kepler K2 Mission schedules two « 75 
day pointings falling within ±10 degrees of the celestial 
equator, and will obtain precise photometric lightcurves for 
thousands of stars covering about 10% of the region of in¬ 
terest. The planned Transiting Exoplanet Survey Satellite 
(TESS) will obtain precise photometric timeseries of 27-day 
duratation for bright FGK-type stars in the entire equato¬ 
rial region. These surveys will allow refinement of a list of 
candidates for GeoSphere comparison stars. 

5.5 Other Problems 

In GSO, thermally-induced geometric changes are domi¬ 
nated by eclipses at equinoxes. We find such orbital changes 

^ For example. The White Sands Test Fa¬ 
cility Remote Hypervelocity Test Laboratory, 

http://www.nasa.gov/centers/wstf/laboratories/hypervelocity 


to be important at the level of 0.001% only for a few minutes 
as GeoSphere enters and exits the Earth’s shadow. 

GeoSphere will be measured only when out of eclipse 
(a > 10.9°) not at opposition. A thin earth crescent is al¬ 
ways visible. Table 1 shows that when separations /3 of the 
Sun and Earth are 20°, the using Lambertian scattering, the 
Earth is 10.8 magnitudes dimmer than the Sun (a flux ra¬ 
tio of 2 X lO'^). While below our target precision, this extra 
light is something to be considered in analyzing GeoSphere 
data. Most likely one will solve for the Sun’s brightness and 
suitable parameterizations of albedo and phase functions of 
GeoSphere and Earth in the analysis. 

For a GSO, GeoSphere observations are best made 
within a few hours of midnight, daily. Irradiance variation 
data show significant power (excluding flares which can be 
identified and removed) on time scales of hours and days, 
that is related in both phase and amplitude to the decadal 
variations. Using hourly irradia nce data from the VIRGO 
instrument (iFrohlich et al.|[l997h . we find that sampling ev¬ 
ery 3 days suffices to reduce high frequency aliasing to 10% 
of the total variance. 

Lastly, care should be taken to determine the polariza¬ 
tion properties of a GeoSphere which is always irradiated 
anisotropically. This can be determined prior to flight and 
using laboratory models, as well as measured directly during 
from telescopes on the ground. 


6 CONCLUSIONS 

We have suggested that an inert scattering sphere be placed 
in a very stable geosynchronous orbit in order to make mea¬ 
surements of the solar flux variations, and measurements 
of the Earth’s changing albedo. There are technical chal¬ 
lenges but no fundamental show-stoppers in making these 
measurements. To achieve a precision better than 0.1% over 
centuries, micrometeoroid-induced surface degradation must 
be mitigated, modeled and/or measured. Such precisions are 
of great interest to both solar and terrestrial climate change 
physics. 

There is perhaps a direct conflict with ITU’s require¬ 
ment for exit strategies on GSO spacecraft. But if the 
community adopts the “stable plane” orbits identified by 
iFriesen et al.l (Il992fl as “graveyard orbits”, the ITU man¬ 
date will not be applicable with a GeoSphere placed in such 
an orbit. 

Finally, one year of ub photometric observations of Geo¬ 
Sphere should reveal phase relationships in spectral irradi¬ 
ance changes above and below 400nm, resolving a debate of 
importance to climate science ([Harder et al.l[200^L 


ACKNOWLEDGMENTS 

We are grateful for discussions with Giuliana de Toma, 
Roberto Gasini, and a conversation with Aaron Rosengren, 
James Garry and Duncan Smith via ResearchGate. The ref¬ 
eree is thanked for very useful suggestions. 


© 2014 RAS, MNRAS 000, [THS] 
























GeoSphere 5 


REFERENCES 

Allen C. W., 1973, Astrophysical Quantities, Athlone 
Press, Univ. London 

Dever J., Banks B., de Groh K., Miller S., 2012, Hand¬ 
book of Environmental Degradation of Materials, Chapt. 
29. Degradation of Spacecraft Materials, Elsevier, Ams¬ 
terdam, 2nd edition 

Frohlich C., 2013, Space Sci. Rev. 176, 237 
Friesen L. J., Jackson, IV A. A., Zook H. A., Kessler D. J., 
1992, J. Geophys. Res. 97, 3845 
Frohlich G., Crommelynck D. A., Wehrli C., Anklin M., 
Dewitte S., Fichot A., Finsterle W., Jimenez A., Chevalier 
A., Roth H., 1997, 175, 267 

Hale G. E., 1908, ApJ 28, 315 

Hall J. C., Henry G. W., Lockwood G. W., 2007, AJ 133, 
2206 

Hall J. C., Henry G. W., Lockwood G. W., Skiff B. A., 
Saar S. H., 2009, AJ 138, 312 
Harder J. W., Fontenla J. M., Pilewskie P., Richard E. C., 
Woods T. N., 2009, Geophys. Res. Lett. 36, 7801 
Harder J. W., Lawrence G. M., Rottman G. J., Woods 
T. N., 2000, in W. L. Barnes (ed.). Earth Observing Sys¬ 
tems V, Vol. 4135 of Society of Photo-Optical Instrumen¬ 
tation Engineers (SPIE) Conference Series, p. 204 
Henry G. W., 1999, PASP 111, 845 
Judge P. G., Lockwood G. W., Radick R. R., Henry G. W., 
Shapiro A. L, Schmutz W., Lindsey C., 2012, A&A 544, 
A88 

Kessler D. J. et ah, 2011, Limiting Future Collision 
Risk to Spacecraft. An assessment of NASA’s meteoroid 
and orbital debris programs, Technical report. National 
Academy of Sciences 

Lockwood G. W., Skiff B. A., Henry G. W., Henry S., 
Radick R. R., Baliunas S. L., Donahue R. A., Soon W., 
2007, ApJS 171, 260 

Lockwood G. W., Thompson D. T., Lutz B. L., Howell 
E. S., 1991, ApJ 368, 287 

Lubin D., Tytler D., Kirkman D., 2012, ApJ 747, L32 
Schwabe H., 1844, Astron. Nachr. 21(495), 233 
Shapiro A. L, Schmutz W., Rozanov E., Schoell M., Haber- 
reiter M., Shapiro A. V., Nyeki S., 2011, A&A 529, A67 
Smirnov V., Semenov A. S., Sokolov V., Konoshenko V., 
Kovalyov I., 2000, Space Debris 2, 1 
Ulivieri C., Circi C., Ortore E., Bunkheila F., Todino F., 
2013, Journal of Guidance Control Dynamics 36, 935 
White O. R. (ed.), 1977, The Solar Output and Its Varia¬ 
tions, Colorado Assoc. Univ. Press, Boulder 
Willson R. C., 2014, Ap&SS 352, 341 
Young A. T., Genet R. M., Boyd L. J., Borucki W. J., Lock- 
wood G. W., Henry G. W., Hall D. S., Smith D. P., Bali- 
umas S. L., Donahue R., Epand D. H., 1991, PASP 103, 
221 


© 2014 RAS, MNRAS 000, [THS] 


